Accidental suppression of Landau damping of the transverse breathing mode in 

elongated Bose-Einstein condensates 

B. Jackson and E. Zaremba 
Department of Physics, Queen's University, Kingston, Ontario K7L 3N6, Canada. 

(February 1, 2008) 



We study transverse radial oscillations of an elongated 
Bose-Einstein condensate using finite temperature simula- 
tions, in the context of a recent experiment at ENS. We 
demonstrate the existence of a mode corresponding to an in- 
phase collective oscillation of both the condensate and ther- 
mal cloud. Excitation of this mode accounts for the very small 
damping rate observed experimentally, and we find excellent 
quantitative agreement between experiment and theory. In 
contrast to other condensate modes, interatomic collisions are 
found to be the dominant damping mechanism in this case. 
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The study of collective excitations provides us with a 
powerful tool for probing the fundamental properties of 
quantum many-body systems. For example, measure- 
ments of low-lying excitations in a Bose-Einstein con- 
densed (BEC) gas of trapped atoms Q helped establish 
the Gross-Pitaevskii equation (GPE) as an excellent de- 
scription of condensate dynamics at low temperatures Q . 
At higher temperatures the situation is complicated by 
the existence of a significant noncondensed fraction. In- 
teractions with this thermal cloud leads to damping and 
frequency shifts of condensate collective modes, as mea- 
sured in several experiments [||-^). A theoretical under- 
standing of these processes is important as a foundation 
for continuing studies in this SXGcl, clS well as providing 
a fascinating opportunity for studying quantum systems 
at finite temperatures. 

One can think of interatomic scattering processes as 
being comprised of mean fields and collisions, where their 
relative importance in describing the condensate dynam- 
ics normally depends upon the rate of collisions compared 
to the characteristic timescale of the oscillation. When 
collisions are relatively frequent, the coupled dynamics 
of the condensate and thermal cloud are most appro- 
priately described by hydrodynamic Landau-Khalatnikov 
two- fluid equations [Q . It is well known that in liquid he- 
lium this leads to hydrodynamic first and second sound 
modes. In the opposite, collisionless regime, mean field 
interactions become more important. With the possi- 
ble exception of Ref . H , experiments have so far resided 
in this regime. Here, coupling between fluctuations in 
the condensate and normal thermal cloud densities leads 
to Landau damping and an associated frequency shift of 
the condensate mode, and has been the subject of sev- 
eral theoretical treatments |^,^|. However, the thermal 
cloud can play a role beyond that of simply providing a 



damping mechanism. It can undergo collective motion of 
its own, so that the response of the system to an external 
probe can give rise to a coupled oscillation of both the 
condensate and thermal cloud. This kind of mean field 
coupling is responsible for some of the behavior 

observed in the experiment of Ref. j| . 

In this paper we study the transverse breathing mode 
of a cigar-shaped condensate at finite temperatures, as 
recently observed in an experiment at ENS As we 
shall see, the coupled dynamics of both the condensate 
and thermal cloud are particularly important in this situ- 
ation, giving rise to somewhat surprising behavior. Such 
considerations are crucial in providing an understand- 
ing of the very small damping rates and frequency shifts 
measured in this experiment. 

Consider a Bose gas with atoms of mass m confined 
in an axisymmetic harmonic trap potential of the form 
V(r) — m[uj^(x 2 + y 2 ) + ujIz 2 ]/2. The ratio between the 
axial and radial trap frequencies, A = lu z /luj_, charac- 
terizes the trap anisotropy. In this work we will apply 
the formalism of |Q , where the dynamics of the system 
are described by a generalized GPE for the condensate 
wavefunction, $(r, t), 

<9<E> / h 2 „ \ 
ih—=l-—V 2 + V + gln c + 2h]-iR\<P, (1) 

and, within a semiclassical Hartree-Fock approximation, 
a Boltzmann equation for the thermal cloud phase-space 
distribution, /(p, r, f), 

% + ~ ■ V/ - VC/ cff ■ V p / = C 12 [f] + C 22 [f}. (2) 
at m 

The condensate and thermal cloud densities are given 
by n c (r,i) = |$(r,i)| 2 and h(r,t) = J(dp/h 3 )f(p,r,t), 
respectively. In the above equations, an important quan- 
tity is the interaction parameter g = Ai:h 2 a/m (where a 
is the s-wave scattering length). It not only enters into 
the effective mean field potential U c s = V + 2g(n c + n), 
but also into the collision integrals C22 (which represents 
binary collisions involving only thermal atoms) and C\% 
(involving a condensate atom in the incoming or outgoing 
channels). The effect of atom exchange in the C12 col- 
lision process appears in (|l|) through the non-Hermitian 
source term R(r,t) = (h/2n c ) f (dp/ h 3 )C 12 [f]. 

First of all, it is instructive to summarize the behavior 
of the system in the T = limit and for temperatures 
above the transition, T > T c . The former corresponds to 
a pure condensate, described by (|l]) with n — and R = 
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0. For an axisymmetric anisotropic trap (A 7^ 1), analysis 
of the GPE in the interaction-dominated Thomas-Fermi 
(TF) limit gives two low-energy m = modes with fre- 
quencies = 2 + 3A 2 /2 ± (9A 4 - 16A 2 + 16) 1 / 2 /2 
M. The high-lying mode lu + will be of interest here, 
which for a cigar-shaped trap with A <C 1 corresponds to 
a breathing-type oscillation in the radial direction, with 
a frequency uj + — > 2lu± as A — > 0. In this limit the radial 
oscillations are decoupled from those in the axial direc- 
tion. Interestingly, for a non-interacting condensate the 
mode frequency is also 2uj± for all A, so that the trans- 
verse mode frequency in the A — > limit is independent 
of interactions. This behavior is similar to that found 
theoretically for a gas in a 2D harmonic trap interact- 
ing with a contact potential [[l3ll , where an underlying 
hidden symmetry gives rise to an undamped breathing 
mode with frequency w = 2luj_, independent of interac- 
tion strength, temperature, or statistics. For the ENS 
experiment, A ~ 6.46 x 10~ 2 , and the high-lying mode 
has the TF frequency u> ~ 2.00052 u±, in good agreement 
with the experimentally measured value. 

Above T c , the thermal cloud evolves according to (|^) 
with n c — and C12 = 0. The mode frequencies can then 
be calculated by taking moments of the kinetic equation 
Ji"4j . In the hydrodynamic limit the transverse mode fre- 
v / T0/3ujj_ IHll for A -> 0, while the 
The ENS experiment 



quency is whd = V ^U/'6uj±_ |14 11 
collisionless result is ojcl — 2u±. 
resides in the near-collisionless regime, where the mode 
has frequency u> ~ luql and is weakly damped through 
interatomic collisions. So, one has the interesting situa- 
tion where both the condensate at T = and the thermal 
cloud at T > T c have mode frequencies close to 2ui±. 

The question then arises as to the behavior for < 
T < T c , where both the condensate and thermal cloud 
are present. Answering this question requires a self- 
consistent solution of the coupled GPE (|l) and Boltz- 
mann equations (^|). We perform this task numerically 
using a split-operator/Monte-Carlo simulation, as de- 
scribed in |16[ ]. As we shall show, a mode exists at finite 
temperatures corresponding to an in-phase motion of the 
two components with a frequency of u> ~ 2luj_, which 
is a direct consequence of the behavior in the T = 
and T > T c limits. This mode, which is only weakly 
damped, corresponds to a breathing oscillation of the 
entire trapped gas. We demonstrate that this mode was 
predominately excited in the ENS experiment [^), lead- 
ing to the anomalously small measured damping rate and 
frequency shift. 

We first discuss results of our simulations for a tem- 
perature of T = 125 nK, where all other parameters are 
chosen to match the experiment ||[l7|. At this temper- 
ature, which should be compared to the experimentally 
measured critical temperature of T c ~ 290 nK, we find 
a condensate fraction of 0.56. We calculate the radial 
moments, (x 2 + y 2 ) = J drn(x 2 + y 2 ), as a function of 
time for both the condensate and thermal clouds. This 
moment projects out the transverse breathing mode, and 



is related to the condensate squared radius measured in 
Ref. §. 

Three distinct excitation schemes are considered. The 
first matches that employed experimentally, where the 
radial trap frequency was abruptly changed, then reset 
to its original value after a time r. This procedure can 
be represented by u' ± (t) = uj±{1 + a[Q(t) — 9(i — r)]}, 
where O(x) is the Heaviside step function. In Fig. 1(a) we 
show results for a = 0.26 and ui±t = 0.172. The figure 
clearly shows that both components respond by oscillat- 
ing with approximately equal amplitudes and phase, at 
a frequency very close to 2u±. The in-phase oscillation 
of the condensate and thermal cloud continues over the 
entire timescale of the simulation, with both amplitudes 
decaying at a very similar rate. Importantly, this damp- 
ing is very weak, much smaller than one would expect for 
Landau damping which is usually the dominant damping 
mechanism in this regime. 
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FIG. 1. Time-dependent moment (x + y ) for the con- 
densate (narrow line) and thermal cloud (bold line), divided 
by the corresponding values at t = 0. Data is for a temper- 
ature of T — 125 nK and is the result of exciting the system 
using the tophat perturbation employed experimentally Q. 
Simulations shown in (a) include collisions, while those in (b) 
exclude all collisions (C12 = C22 ~ 0). 

To understand this behavior one should note that cal- 
culations of Landau damping [^,^) conventionally assume 
that the noncondensate is always in thermal equilibrium. 
Fluctuations of the condensate mean field induce a re- 
sponse in the thermal cloud density, where the dissipative 
and reactive components in turn give rise to damping and 
frequency shifts of the condensate mode. However, this 
picture is clearly inadequate when the thermal cloud is 
itself in oscillation. In this case, mean field coupling may 
lead to the thermal cloud driving the condensate at its 
own resonant frequency, which can be significant when 
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the natural frequencies of the two components are nearly 
degenerate. As shown in [fj"l|| , a similar near-resonant 
driving of the condensate is in fact responsible for the 
abrupt change in mode frequency observed at high tem- 
peratures in Ref. ||. The present case can be seen as 
an extreme example of this behavior, as the natural fre- 
quencies of the condensate and thermal cloud are both 
very close to 2wj_. Mean field interactions then lead to 
phase-locking of the oscillations of the two components. 

We test this interpretation by separately exciting the 
two components in our simulations. This is readily 
achieved by imposing velocity fields at t = of the form 
v ex xi + yj 1 1 ll| on either the condensate or the ther- 
mal cloud. Fig. 2(a) shows the result of exciting the 
condensate only. The condensate oscillation is initially 
damped quite strongly, while the oscillation in the ther- 
mal cloud slowly builds up in time. Eventually, the con- 
densate and thermal cloud settle down to a small ampli- 
tude in-phase mode. The initial condensate relaxation 
is essentially Landau damping due to mean field interac- 
tion with the stationary thermal cloud, with a damping 
rate similar to that found for the m = 2 mode discussed 
later. This excitation scheme could be implemented ex- 
perimentally by optically imprinting a phase onto the 
condensate, and would be useful in verifying the physical 
picture presented here. 
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FIG. 2. Results of initially exciting only the (a) conden- 
sate, and (b) thermal cloud by imposing velocity fields on 
each, for a temperature T = 125 nK. In (a) the condensate 
moment is represented by the narrow line with the thermal 
cloud in bold, while in (b) the line styles are interchanged for 
clarity. 

In Fig. 2(b) we show the result of exciting the thermal 
cloud only. Again, we see the relaxation of one compo- 
nent (thermal cloud) matched by excitation of the other 
(condensate), with both settling into an in-phase oscilla- 



tion at later times. The larger amplitude of the in-phase 
mode in this case is a consequence of the much greater 
spatial extent of the thermal cloud. On imposing the 
same velocity field, the energy increase is proportional to 
the equilibrium value of (a; 2 + y 2 ), which is several times 
larger for the thermal cloud than for the condensate at 
the temperature considered. During the subsequent evo- 
lution, the initial energy of either component is redis- 
tributed until the two components oscillate together in 
phase. Since the energy of this in-phase mode is mostly 
stored in the thermal cloud, one can see that the final am- 
plitude should indeed be larger when the thermal cloud 
is excited initially. 

To compare our results directly to experiment, we fit 
the condensate data for our first excitation scheme [as 
displayed in Fig. 1(a)] to an exponentially decaying sinu- 
soid A cos(ujt + (p)e~ rt + C. Our results for the frequency 
u>, and damping rate T Jl9| are plotted in Fig. 3 for a 
range of temperatures, together with the corresponding 
experimental data. We find very good agreement be- 
tween experiment and theory, for both the frequency and 
damping rate [Q. It is interesting to note that, along 
with the very small damping, the frequency is virtually 
independent of temperature in both simulations and ex- 
periment. This is in marked contrast to experimental 
|5|-|5| and theoretical results flp||llll for other modes, and 
is a unique feature of this breathing mode excitation. We 
note that fits to the thermal cloud data give results very 
similar to those in Fig. 3. 

So far we have not addressed the origins of the damp- 
ing. To provide insight into this question, we have per- 
formed a simulation at T = 125 nK where collisions are 
artificially switched off. Thus, the only coupling between 
the two components is provided by the mean field terms. 
The result is plotted in Fig. 1(b). The most striking 
feature of this plot, when compared to Fig. 1(a), is the 
almost complete absence of damping, either in the con- 
densate or the thermal cloud. Specifically, a fit to the 
condensate moment yields a damping rate an order of 
magnitude smaller than the collisional result. Further- 
more, the separate inclusion of cither C12 or C22 colli- 
sions yields values intermediate between the collisionless 
and full collisional results. Finally, artificially increasing 
the collision cross-section by a factor of two leads to an 
increase in the damping rate by a corresponding factor. 
These simulations, together with the excellent agreement 
found with experiment, provide strong evidence that col- 
lisions play a crucial role in determining the observed 
damping. 

We should emphasize that the supression of Landau 
damping discussed here is specific to the m = trans- 
verse breathing mode in a highly prolate harmonic trap, 
and is a consequence of the accidental degeneracy be- 
tween the condensate and thermal cloud frequencies. 
This condition is no longer satisfied, for example, by the 
lowest m — 2 mode, since the condensate now has a 
frequency of approximately \/2u)±_ while the thermal 
cloud frequency remains close to 2u>± . A simulation 
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at T = 125 nK for this mode indeed shows a damping rate 
of r ~ 50 s^ 1 , an order of magnitude larger than that of 
the m = mode and in agreement with the qualitative 
result mentioned in Q. Moreover, we would also expect 
larger damping rates for the m — mode as A — > 1, since 
the condensate TF frequency approaches \fhuu±_ in this 
limit. Simulations at A = 0.57 and A = 0.75 confirm that 
lifting the degeneracy between the condensate and ther- 
mal cloud modes leads to an increase in damping with 
decreasing trap anisotropy. 
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FIG. 3. Frequency (a) and damping rate (b) of the con- 
densate breathing mode. Our results (open circles) are com- 
pared to the experimental data (solid circles), where the sim- 
ulation parameters and excitation scheme are chosen to re- 
produce the experimental conditions || . 

In summary, we have used hnitc-tcmpcrature simula- 
tions to study the transverse breathing mode in elon- 
gated Bose-Einstein condensates, as realized in a recent 
experiment We show that the very small damping 
rate and frequency shift observed experimentally is a re- 
sult of a degeneracy between the condensate and ther- 
mal cloud oscillation frequencies, which results in a mode 
comprised of an in-phase oscillation of both components. 
The weak damping arises through interatomic binary 
collisions, rather than mean field interactions (Landau 
damping) as is normally the case. A comparison of sim- 
ulation results to experiment shows excellent agreement, 
confirming the validity of our theoretical approach and 
underlining the importance of treating the full noncon- 
densate dynamics, including collisions, when accounting 
for experimental behavior. 
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In particular, we consider 87 Rb atoms confined with a ra- 
dial trap frequency of ui±_ — 2n x 182.6 Hz. Experimental 
use of evaporative cooling meant that the total number of 
atoms decreased with temperature. We account for this 
variation in our simulations. 

Note that in the limit lj±t <C 1 the tophat excitation 
employed experimentally corresponds to a delta function 
impulsive perturbation. In linear response this is equiva- 
lent to imposing identical velocity fields of the form given 
on both the condensate and thermal cloud. 
Due to small departures of the condensate data from an 
exponentially-decaying form, our determination of the 
damping rate has some uncertainty, which we estimate 
to be about 15 %. 

Experimentally a condensate amplitude revival was ob- 
served at low temperatures and after strong excitations. 
We see no evidence of this behavior in simulations that 
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extend over periods twice as long as those in Fig. 1. Pos- 
sible sources for this effect are discussed in H. 
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